ON A BROWNIAN EXCURSION LAW, I: 
CONVOLUTION REPRESENTATIONS 



Michael Schroder 

Abstract: This paper studies Brownian motion subject to the occurrence 
of a minimal length excursion below a given excursion level. The law of this 
process is determined. The characterization is explicit and shows by a layer 
construction how the law is built up over time in terms of the laws of sums 
of a given set of independent random variables. 
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1. Introduction 



In this paper we identify a general Parisian-type excursion law and determine its structure 
in the Brownian case. The idea is that we do not just ask when a process passes a given 
level, but instead ask for how long it will thereafter stay on the same side of the level. 
These excursions may have any length and may occur at any time. Having settled on a 
minimal length for them, the question to ask is if a minimal-length excursion will occur 
during a given time span. This then yields a second source of randomness for the problem, 
and the object to be studied is its joint law with that of the given process; see Section 2 
for how to make all this precise. 

The principal difficulty is to obtain an explicit description of this joint law when the process 
to start with is given explicitly. The main contribution of the paper is to give such an 
explicit description when this process is Brownian motion; details are in Section 3. 

Historically, a rigorous study of these questions seems to originate with [13]. Extensions 
^ \ and further applications are developed for example in [6] and the references there. Appli- 
cations to insurance are emerging; see for example [5]. 

Regarding the structural understanding of the joint law, its Laplace transform with respect 
to time has been determined in the Brownian case; see [13, Appendix] with modifications in 
[11]. As made precise in Section 5, the Laplace transform permits one to separate the above 
two sources of randomness of the problem. This is effected in [13, Appendix] by recourse 
to the Brownian meander and the Azema martingale; the key steps of the argument are 
recalled in Section 6. The result is quotients of higher transcendental functions as Laplace 
transform of the joint law, with the denominators corresponding to the law's excursion 
source of stochasticity and the numerators to the law's process source of stochasticity. 
This generalizes situations as encountered for example with Laplace transforms of first 
passage times of Brownian motion, which lead to laws in terms of theta functions. 

Mathematically, the contribution of the paper is to provide analogous such representations 
to the present context. The basis for this is furnished by our finding in Appendix B that the 
above denominators satisfy a functional equation, and as its main technical contribution 
the paper turns this into a way for the analytic inversion of the corresponding quotients. 
With more details developed in Section 4, this then shows how the effect of the excursion 
source on the law is built up over time - by adding at any integer time n + 1 a new layer 
of n independent copies of a fixed explicit random variable. This picture is made rigorous 
in Section 9 yielding the explicit expressions for the joint law of Section 3. 
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To conclude, Section 10 indicates how our approach is complemented in the companion 
[12] to the present paper by constructive techniques for convolution representations. 

2. Basic notions and facts: the general process excursion law 

We address in the paper the occurrence of excursions of a minimal length duration of 
a given stochastic process below a given excursion level which take place within a given 
period of time, and ask to characterize the law of the new stochastic process thus obtained. 
We follow [11] to make this precise in two stages. A general framework is thus established 
in the present section which furnishes the basis for the Brownian case to be developed 
subsequently in the paper. 

2.1 Basic setting: Let X be any real- valued stochastic process on the complete proba- 
bility space (O, J 7 , Q) whose time set is the nonnegative reals [0, oo), so that X = (X u ) u > , 
and let it be adapted to the filtration F = {^u)u>Q on. this space which satisfies the usual 
conditions. 

As the primitives of the problem fix any time t, specify any real a, which is to take the 
role of the excursion level of the problem, and specify any real D > 0, which is to take the 
role of the minimal excursion duration. 

Denote by T t>a := T(X) tjCL the first passage time after time t to the level a of X given by 
T(X) tja = inf{s > t\X s G [a, oo)}; for simplicity assume that T tjCL is a stopping time, and 
if t = we drop reference to the time t from the notation. 

2.2 Achievement time: We first recall how the excursion situation under considera- 
tion is encoded by the achievement time H a t := H(X) ajt , a random variable to be defined 
presently with the property that for any real T > t we have H t a < T iff there is in [t, T] a 
subinterval of length at least D on which X only takes values below the level a. Morally, 
H a: t so is the upper bound of the first of these subintervals occurring after time t. For its 
explicit construction, which is to follow, for any real u > introduce a ajU = a(X) a>u to 
denote the smallest upper bound of subintervals as above of [u, oo), namely: 

a a>u = inf{s G [u+D, oo) | A| (s _ D)S) < a} . u G [0, oo) . 

We will need to distinguish the two cases when X t is equal to or above the level a, and 
when it is below this level. 

If X t > a, the situation to be referred to as the first case, the question arises if X will 
become and remain smaller than a for all points in time of an interval of length at least D 
contained in [t, oo]. We represent this first case by setting: 

H a ,t = cr a , t if X t > a . 

If X t < a, the situation to be referred to as the second case, a new idea is needed when 
X has been staying below a for some connected period of time already. This excursion 
needs to continue below a for a period of length 5t G (0, D) only to reach the minimum 
duration D stipulated. Two subcases result which we formalize in terms of the location of 
t+5 t relative to the first passage time T^ a = T(X) t ^ a . Firstly, if T^ a > t+5 t , the situation 
to be referred to as Subcase 1 , then the process X continues to stay below a during the 
period of time from time t until time t+5t a fortiori, whence the representation: 

H a , t = t+5 t if X t < a and T M > t+5 t . 
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Secondly, if Tt, a < t+dt, the situation to be referred to as Subcase 2, the level a is passed 
by X earlier than time t+St, and the clock for reaching the minimum length D is restarted 
at time T^ a . This then puts us into the first case, albeit as of time T tja instead of time t, 
whence the representation: 

H a , t = <r a ,T t , a if X t < a and T tja < t+S t . 

2.3 Excursion density: The excursion law to be studied is then by definition the 
representing measure of the functional given by: 

>->• E [l{H a , t <T} H X T) I t] , 

for any Q-measurable function on the reals R; assuming absolute continuity with respect 
to Lebesgue measure in addition the concept to be studied is hence the excursion density 
h a ,t '■= h(X) at t which by definition is the conditional density on [t, oo) x R characterized by: 

E[l {HaX t)O(Xt) -V = J* R <j>(x) h a j(T, x) dx , 

for any Q- measurable function (p on the reals R; here notice x h-» h ajt (T, x) = for T < 5 t 
extending the definition of St in Section 2.2 to the first case there by the convention St '■= D. 

2.4 Normalizations: We develop a normalized form of the excursion law. The normal- 
ization is based on an emulation of the restarting at a fixed stopping of Markov processes. 
The construction starts from postulating the existence of a process X* which satisfies 

X* u l ^X t+u -X t , ue [0,oo), 

as well as X*(0) = 0, and which is independent of J^; examples for processes X for which 
these normalized processes exist thus include Brownian motion, Levy processes which creep 
across levels, or continuous Markov processes which are homogeneous in their time and 
their space variable. The construction of Section 2.2 applied to X* then yields for any real 
b the normalized achievement time H£ := H(X*)^q which satisfies 

H a,t = t + H*_ Xt , 

as well as 

E[l {Hat < u+t} cf>(X u+ t)\jr t ] = E[l {H: _ x < u} <f>(X t +X:)} , 

for any measurable map (f> on R, treating the time-t value X t of X as a real number. 
Granting absolute continuity with respect to Lebesgue measure in addition as in Section 2.3 
above, densities hi on [0, oo) x R are hence defined by: 

E[l {H ;<u } r(K)] = f R cf>*(y)K(u,y)dy, 

for any Q-measurable function (jf on R. By construction these normalized excursion den- 
sities hi satisfy hl([0,So) x R) = {0} with Sq the minimal excursion duration remaining 
as in Section 2.3 above, and are related to the excursion densities h a ,t of Section 2.3 by: 

h a , t (u+t, y+X t ) = K_ Xt (u, y) , 

for any real u > and y. Granting existence, the study of the former law is hence reduced 
to the latter one. Translating the first and second case of Section 2.2 this explicitly asks 
to characterize hi, for any real b, in the following Cases I and II respectively: 

(I) We have b < and ask if there will be a length D time interval Id on which X* | Id < b. 
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(II) We have b > and ask if X* | [ (5o ) < b where the positive real 5q < D is the Section 2.2 
minimal excursion length remaining. 

Here Case II further decomposes into two subcases according to the relative position of Sq 
to T 6 *, the first passage time of X* to the level 6, as follows. If moreover T 6 * > So, we are 
in subcase (II-l), otherwise we are in subcase (II-2). 

The paper initiates a study of the structure of hi along these lines when X is Brownian 
motion; our main results here are to be described in the next Section 3 as a first step. 

3. Statement and discussion of main results: 
convolution representations of the excursion law in the Brownian case 

In this section we formulate and discuss our main results about the normalized excursion 
law hi of Section 2.4 in the Brownian case. We delineate in Sections 3.2 and 3.3 how, 
with hindsight, it emerges from a core relationship pinned down in Section 4 by taking 
into account the additional interrelations entailed by the Case I and Case II situations. 
Our discussion is in terms of convolutions of laws, with representing functions listed in 
Section 3.1 to follow. 

3.1 Set up and functions: We modify the Section 2.1 set up as follows. First we let the 
probability space (O, J 7 , Q) there be equipped with the standard filtration F = | u G 
[0, oo)) of the Brownian motion W on it. Then let X = W and, consequently, let X* = W* 
in Section 2.4 with W* now being the standard Brownian motion independent of time-t 
information J^~ t that is obtained by time-t restarting of W. 

It is convolution representations for hi which thus result. They are in terms of three main 
classes of functions. Referring to Appendix A for more detail about them, these are as 
follows. Firstly, the functions u n on (0, oo) given for any integer n > 1 by the n-fold 
convolutions on (0, oo): 

v n = where u(u) = (2/^)^/(2^ + 1), u E [0, oo) , 

and by vq = v^°> being set equal to the Dirac delta function at 0; for n > 1 note v n = 
law (nV) with nV the n-fold sum of independent copies of a random variable V with 
law(y) = v. Secondly, for any real a and c, the functions p a , c on (0, oo) given by: 

_ / ^DX\a+c\ C < 0, 

Pa,c — \ /-=- 

[VDx\a\*9c c > 0; 

they are convolutions on (0, oo) in terms of two additional functions on (0, oo) that depend 
on the complex parameter a with |arg(a)| < j and for any real u > are given by: 
X«(«) = (l/\fwu ) exp ( — (a/2) 2 /u) respectively: 

Sa (u) = exp (- ^){ ft(M )exp(- ^£^) + ^Erfc 

Thirdly, for any excursion level b the function (3 = (3^ on the real line R given by: 

(5{y) = {b-y)/VD. 

3.2 Case I results: We first assume the situation of Case I to hold as it is described 
in Section 2.4; we are thus asking about the occurrence of an excursion of duration at 
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least D below the level b < at some point in the future. The structure of the normalized 
excursion law hi here is governed by the sign of (3 on the 'state space' by way of the 
functions p a , c - Our precise result is as follows. 

Theorem 3.1: In the above Case I setting we have for any reals u > and y the convo- 
lution sum representation: 

y) = w (""•«»> * "* (n_1) ) ( (1/2) ^ " "))• 

where the sum is over all integers n > 1 satisfying n < u/D, and where b* = b/y/D. 

With a proof in Section 9.1 our results make explicit how the law is built up over time 
by the addition of a new layer at integer points in time. Anticipating the discussion in 
Sections 5 and 6 this is a principal feature. It gives expression to the effect of X\ = 
as the first of two sources of stochasticity that determine the structure of the excursion 
law. The effect of the second source of stochasticity X 2 = W*(H£), on the other hand, is 
encoded by the structure of the layers to be added. Here it gives rise to a fixed density 
P = Pb*,p(y) which, at each time n, is convolved with v n -i, the (n — l)-fold convolution 
of the density v originating with the above first source of stochasticity. The structure of 
the function p depends on the state variable y, and gives expression to the position of y 
relative to the excursion level b by way of the function f3. 

3.3 Case II results: We next address the structure of the normalized excursion law 
hi in the situation of Case II of Section 2.4; recall that this setting is characteristic for 
the problem by asking about the excursion of Brownian motion W* to continue below the 
level b > for a period of length at least 5 := 5 > smaller than D. We here obtain hi 
as a four term sum: 

K = Sj=i J2k=i K,j,k > 

where the single summands are functions of the time variable u > and the state variable 
y as follows. Pertaining to the situations where the first source of stochasticity X\ = H£ 
(see the discussion of Section 3.2) is in Case II- 1, 

Q*(T b >5)Q*(T b <5) l i5tQo) (u) r oo_ v ^_^ 2 (p-x-y)^ 
lM>,?/) ~ D y/2^S) Jo H 2D 2(u-S) r*' 

K,i,2( u ^y) = Q*(T b > 5) i(s,oo)(u)<Pb,u(y); 

here cp b ,u is a translate of the transition density from to y < b of Brownian motion killed 
at the first hitting time of b and living on (— oo, b) as follows: 

where Q*(T b > S) = (l-Erfc)(6/v^) and Q*(T b < 6) = Erfc(6/v / 25), and where the 
integral factor of hi 1 1 is expressible in terms of Erfc and its derivative as well; see the 
representation for the function h b ^ in [11, p. 4]. 

The remaining two summands have been known at the level of their Laplace transforms. 
By inverting these transforms in Sections 9.2 and 9.3 we are now able to describe in the 
next two results how their structure is built up over time. This once more proceeds by 
expressing them as sums of convolutions in terms of Section 3.1 functions. Our first result 
pertains to the summand h* h21 as follows. 
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Theorem 3.2: In the above Case II situation, with excursion level b > 0, we have for 
any reals u > and y the convolution sum representation: 

I* / \ 1 r SA(u-Dn) I \ *(n-l)\ (U-Dn-W\ 

n ^ ' 

where the sum is over all integers n > 1 satisfying n < u/D. 

With a proof in Section 9.2 this theorem refers to a situation where the first stochasticity 
source Xi = H£ is now considered on the set of events where is expected to be restarted. 
This is given expression to by integration with respect to [i b , the law of the first passage 
time of Brownian motion to the level b given by: 

Li b (dw) = i) bV z(w)dw; 

see Section A.l. Conditional on that influence, the excursion law is built up at positive 
integer time points as follows. At each time n add the new term obtained by convolving 
z/ n _i, the (n — l)-fold convolution of the density v, with a fixed density p b (-,y). Here the 
function p b on (0, oo) x R is explicitly given by the following integral: 

Pb(r, y) = 2 J R (p b:S (x) N 0iDT (\x-y\/V2) dx , 

where N 0jV (^) = (l/y/2nv) Jj-oo^] exp (— x 2 /(2v)) dx denotes the normal distribution with 
mean equal to and variance equal to v. The functions p b originate with the restriction of 
the second source of stochasticity X 2 = W*(H^) to the set of all events where a restarting 
of HI will not happen. 

The situation for the final fourth summand h\ 22 differs from the above set up in regard to 
the stochasticity source X 2 = W*(H^): here it is considered on the set of all events where 
a restarting of its argument is to take place. Conditional on the restarting of Xi as 
above we therefore are in the situation of Section 3.2, and hence seek to characterize the 
occurrence of future length D excursions below the level 0. This situation is encoded by 
the Section 3.1 functions Po,p( y ) and our precise result is as follows. 

Theorem 3.3: In the above Case II situation, with excursion level b > 0, we have for 
any reals u > and y the convolution sum representation: 

K,2,2(u,y) 

Q*(T b < 5) ^ (-l)™" 1 r SA(u-Dn) *(n-l)\f u - Dn - w \ 

n v ' 

where the sum is over all integers n > 1 satisfying n < u/D. 

This result is proved Section 9.3, as the final step of our argument which starts in Section 4. 

4. The key result for Laplace inversion 

At the heart of the structure of the normalized excursion law hi as expressed by the results 
of Section 3 is the reconstruction of the law of Hq, the level-0 normalized achievement time. 
This reconstruction is to be established in this section, and asserts Hq morally to originate 
as (the weak limit of) an infinite sum of independent copies of a fixed random variable, to 
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be denoted by V, in which only finitely many summands contribute when looked at in a 
pointwise sense. Here V is pinned down in terms of its distribution by: 

law(F) = v where v(u) = (2/y/ir) Vu/(2u + l) , u G [0, oo); 

our construction, however, necessitates two auxiliary random variable summands as follows. 

Theorem 4.1: For any real u > 0, we have the equality of laws of pairwise independent 
random variables: 

lw(A R +H*)(u) = E x < n<2n \^j2 l ^( B R + (n-l)V)(u- \n) , 

where A R and B R are random variables such that on some complex half-plane {Rez > zq} 
with zq > we have: 

E[exp(-zA R )] = R(z) and E[exp(-zB R )} = R{z)/s/z , 

for an analytic function R which for some a > \ satisfies R{z) = 0(\z\~ a ) as \z\ — > oo. 

The starting point of this result is provided by joining the equality 



E[ex.p(-zH*)] = l/^(v / 2~Di) , Re (z) > , 

anticipated from Section 6.4 in terms of the function ty(w) = /(o,oo) ^exp (— x 2 /2+wx) dx 
of Appendix B, with the symmetry property of Hq expressed by the functional equation 



\&(w) = V(-w) + v / 27rwexp(i«) 2 ' 

for any complex w, furnished by the key identity of Appendix B. This enables a represen- 
tation of the expectation by resolution as a geometric series: 

= R(z)J^L = 4il -till exp ( - i(„+l),)iV lW », 

where f(z) = (1/V2irz)e~ z/2 and p(z) = -f(z)^(-^/z), where Nx(z) = ty(-y/z)/y/z 
from Appendix A, and granting that for some ao > the first two functions are smaller than 
1 in absolute value when Re (z) > ao. Granting too that Laplace inversion of this sum can 
then be effected term by term, the assertion of Theorem 4.1 follows noting i? 7_1 (A^{ 1 ) = v n 
from Proposition A.l. With the analytical facts implied by the Appendix B leading term 
expansion of omitting further detail the proof of Theorem 4.1 is complete. 

Remark 4.2: The proof of Theorem 4.1 can be seen as being effected by transfer to the 
framework of the Ito theory of excursions of Brownian motion. Denoting the objects of this 
theory by e, this perspective then in particular affords an interpretation of the function p 
at the heart of our inversion in terms of volumes of the Ito measure n as follows: 

p(V2DX) = (V2/2)— A>0 , 

n({e > containing a A-marked point}) 

setting <£>£>(£) = l{t>r>} n({e > of lifetime(e) > £}), for any real t > 0. 

5. Laplace transforms of the excursion law: principal results 

In this section we identify a principal structure of the excursion law enforced when working 
with processes which are restartable twice in the sense of Section 2.4. The thrust is to seek 
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giving expression to the effects of the two sources of stochasticity featured in Section 3 by 
admitting the use of transform methods in time direction, and to identify a framework of 
sufficient conditions for rendering this rigorous; this framework is then to be verified in 
particular when working with Brownian motion in the further development of the paper. 

We therefore concentrate now on the Section 2.4 normalized excursion laws hi associated 
with Markov processes X restartable-at-a-stopping-time twice (with the processes that 
result denoted by X* and X** = (X*)* respectively), and, for simplicity, assume absolute 
continuity with respect to Lebesgue measure of all random variables under consideration 
to hold as well. For the study of the existence of hi we adopt the Laplace transform with 
respect to time, 2?; see Appendix A for pertinent notation and concepts. 

5.1 Basic Laplace transform: Our starting point is furnished by the representation 
of h* h in terms of the transition densities of X as follows: 

K( u i V) = E[1 { h;< u }Px{u-H^ X* H *,y)], 

for any reals u > and y, where pj(r, £, dx) denotes the density of X conditional on time- 
r information subject to £ = X t . To obtain this representation write in Section 2.3 the 
expectation defining h\ as an iterated one by taking it conditional on time-i^ information 
first, and then express these conditional expectations by time-Hj* restarting of X* in terms 
of the transition densities px- 

Two sources of stochasticity are thus identified to determine the structure of hi, namely 
H£ and X*(H£). They can be separated on taking Laplace transforms with respect to 
time once their independence is granted, and the following result is seen to hold then. 

Proposition 5.1: Let the Markov process X be twice restartable and assume H^ and 
X*(H£) to be independent. Then we have for any real y the Laplace transform identity: 

^(hl(-,y))(z) = E[exp(-zH* b )]E[^( Px (-,X* HV y))(z)], 

for any complex z with Re (z) > zq, in the sense of measurable functions with both sides 
either oo or finite. 

As a next step we distinguish the effects of the Section 2.3 Cases I and II on H^. In (II-l) 
observe H£ = Sq. In (I) and (II-2) a reduction occurs to a time-0 and excursion-level-0 
case I situation by way of the decomposition of H£ into independent random variables: 

H b = 1 b + H 0,0 

which gives expression to the very construction of H£ in these follows: the con- 

struction refers to the process X** obtained by restarting of X* at time T & *, the first 
passage time of X* to the level b, and starting from T b * it proceeds by measuring Hq * , the 
achievement time associated with X** by the Section 2.2 construction in the case I situation 
where a = and t = there. Observe that T b * is equal in law to the first passage time 
of X to the level 6, and Hq* is equal in law to Hq, the Section 2.4 normalized achievement 
time for X* with 6 = there. Effects thereof are explained in Sections 5.2 and 5.3. 

5.2 Case I specializations: First consider the Case I situation as left in Section 5.1. 
Using the independence of the summands in the decomposition of H£ derived there, the 
Laplace transform of Proposition 5.1 is here checked to become a product of three Laplace 
transforms as follows. 
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Proposition 5.2: Let the Markov process X be twice restartable. Assume a Case I 
situation in which the random variables H b and X*(H b ) are independent and the Laplace 
transforms J^ ? (hl(-,y)) for any real y are well-defined on the half-plane {Re (z) > zq} 
within the right-hand complex half-plane. Then, for any complex z with Re (z) > zq, 

/(/'/>• V))(z) = E[e- zH ° ] E[e-"S }E[^(p x (; b+X* H *,y))(z)} . 

5.3 Case II specializations: Concentrate on the Case II situation as left in Sec- 
tion 5.1. Referring to Section 2.3 we hence ask about the excursion of X* to continue 
below the level b > for a period of length at least 5 := 5q which is positive and smaller 
than D: 6 G (0,D). The principal structure of the Laplace transform of the normalized 
excursion law hi here mirrors an additional characteristic feature: the relative position of 
5 to T 6 *, the first passage time of X* to the level b. This affords a decomposition of the 
sample space into the disjoint subsets A\ = {T b > 5} and A 2 = {T b * < 5}, corresponding 
to the Subcases II- 1 and II-2 of Section 2.4 respectively. In Proposition 5.1 factors of the 
Laplace transform thus decompose into two term sums. A four term decomposition of 
the Laplace transform results whose summands are Laplace transforms as well. We hence 
obtain four functions hi ■ k on [0, oo) x R, indexed by j, k in {1, 2}, by way of the Laplace 
transform identities: 

&(ht ij , k (;y))(z) = E[l Ai e W (-zH* b )]E[l Ak ^(p x (;X* H; ,y))(z)], 
for any complex z with Re (z) > zq. The precise decomposition result is as follows. 

Proposition 5.3: Let the Markov process X be twice restartable. Assume a Case II 
situation in which the random variables H b and X*(H£) are independent and the Laplace 
transforms J5?(h%(-,y)) for any real y are well-defined on the half-plane {Re (z) > zq} 
within the right-hand complex half-plane. Then the functions hi ■ k are well-defined and 
afford the four-term decomposition: 

K = Sj=i J2k=i K,j,k ■ 

Explicitly we here have: 

K,i,i( u >y) = Q( T b > l{u>5}E[l {M (x*) s <b}Px(u-5,XI,y)] , 

K,i,2(u,y) = Q(T b * > 6)l {u>s} E[l {T .< sy px(u-6,b+X^ 5 ,y)] , 

for any real u > and y, using the running maximum M(X*)$ = max{X* | r G [0,5]}. 
The remaining two functions are for any real y characterized by the three-factor Laplace 
transform identities on {Re (z) > zq} as follows: 

^W,2,i(-,y))(*) = E V~ ZHI } E[l {T -<6}e-* r i } E[l {M{x * )s<b} ^(p x (;X* 5 ,y))(z)], 

^(h* w {-, y )){z) = E[e~ zH ° ] E [1{t*<<$} e~ zT *' ] E[l [T; < 5} ^{p x (;b+X* H;) ,y))(z)]. 

This result is essentially extracted from [13, Appendix Section 8.3] and ultimately makes 
explicit the effects of the Section 5.1 restarting at T b * decomposition of H b once more. 
Seen in conjunction with Proposition 5.2, rendering H b explicit thus reduces to making 
explicit the law of Hq and X*(Xq). We look at this in the Brownian case. 

6. Review: Laplace transforms in the Brownian case 

The general structure of the Laplace transforms of hi has been identified in Section 5. We 
now revert to the Brownian setting of Section 3, and as a second step in the proof of the 
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Section 3 results recall how to make the Section 5.2 and 5.3 Laplace transforms of h\ 
explicit when X = W* is Brownian motion there. The argument follows [13, Appendix]; 
it is based on the structure theory of the Brownian meander and the Azema martingale as 
developed in [1] and [2], with another exposition in [14, Section 12.3]. 

6.1 Brownian motion at normalized achievement time: As a first step in ob- 
taining the Laplace transforms this section reports how the assumptions of the Section 4 
results are satisfied in the Brownian case as follows. 

Fact 6.1: For any real b, the random variables HI and W*(H£) are independent, and 
for the law of the absolute value of the latter we have on R: 

\Wff*\(dx) = l(o,oo ) (x) x exp ( - \x 2 ) dx . 

Simplifying notation B = W* we recall how this is based on the Brownian meander which 
for any real T > is the process mr given by 

tot(u) = \W*(g T + u(t-g T ))\/y/ T -9T, ue [0,1]; 

Here gr = sup{s | s < T and B(s) = 0} is the last time before T where B is 0. The idea 
is to study m T with T = H£. Since ,9~ gT C J^ T by [10, XII (3.2) Lemma, p. 464], it is 
sufficient to show that mr is independent of ^~ 9T and sgn(£? T ), and has the law indicated 
independently of u. Since gr = Tgi by Brownian scaling, a further reduction occurs to the 
case T = 1 where these independence properties are taken care of by the following result. 

Fact 6.2: With respect to the measure Q* = (\Bi\/ E[\Bi\])Q , the process mi is a dimen- 
sion 3 Bessel process BES(3) which is independent of ,W gi . 

In [2, §Theoreme, p. 293] or [14, Section 12.3.2] this is proved as a consequence of a 
generalized Girsanov argument; the idea is to thus construct a Brownian motion (3 such 
that mi satisfies mi (it) = /?(«) + / (l/mi(s))ds as it characterizes Bessel processes of 
dimension 3 started at at time 0. 

As indicated in [2, p. 294] or [14, Step 3, p. 45], the independence properties of Fact 6.2 
entail a relation going back at least to [7] as follows. 

Fact 6.3: We have law(mi)(u,a;) = v^T^" 1 law(BES(3))(l, x) for any u in [0, 1] and 
any real x > 0, where law(BES(3))(l,x) = 2~ 1 / 2 r(3/2)- V exp (-x 2 /2) . 

6.2 Case I Laplace transforms: Keeping the concepts and the notation of Sec- 
tion 6.1, we look at the consequences of the results there for the Laplace transforms of the 
normalized excursion law h£ in Case I, and indicate how a reduction occurs to the Laplace 
transform of the Section 2.4 normalized achievement time Hq. Indeed, as a consequence 
of Fact 6.1 the independence assumptions for Proposition 5.2 are satisfied. The law of 
Brownian motion at time Hq as it enters into the Laplace transform representation there 
is explicitly known by Fact 6.1 as well. Summarizing the computations of [11, Sections 8 
and 9] following [13, Appendix Sections 8.3.1 and 8.3.2], we thus have for any complex z 
with Re (z) > the representation: 

^(h* b (-,y))(z) = E*[ey. V {-zH*)]R I {y,V2D~z), 
in the sense of measurable functions, where the factor 
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Rj(y,V^) = E*[e W (-zT b *)]E*[^( X \x 2 -y\)(2z)] 
= exp 

with b* = b/y/D < 0, is well-defined and finite on {Re (z) > 0}. Here recall the functions 
Xa from Section A.l, G a from Section A. 3, and f3 from Section 3.1. 

6.3 Case II Laplace transforms: Keeping the concepts and the notation of Sec- 
tion 6.1, we study the Case II situation. Here a reasoning analogous to that of Section 6.2 
yields the Laplace transform representation of h b of Proposition 5.3 in the sense of mea- 
surable functions. A reduction occurs of this law to that of determining two Laplace trans- 
forms. The computations of [11, Section 10] following [13, Appendix Section 8.3.3] yield 
for the the first of these, the Laplace transform of the density /jj 21 , the representation: 

&(h* bt2tl (;y))(z)= f* exp(-zw) E*[exp(-zH^]R 2A (y,V2D^)n b (dw) , 
for any complex z with Re (z) > 0, in the sense of measurable functions. Here the factor 

R 2A ( yi V^z)^E*[l {T:>s} ^( XlX2 _ yl )(2z)] 

= VDf n ^( Xlx _ yl/ ^)(2Dz) Vb , s (x) dx 

is well-defined and finite on {Re (z) > 0}, and the remaining concepts are as follows: from 
Section A.l we have the density /i b of T b , the first passage time of Brownian motion W* 
to the level 6, as well as the functions Xa, while the function (p bj g is from Section 3.3. 

For the second Laplace transform in question, the one of the density ^ 21 , we have for 
any complex z with Re (z) > the representation: 

&(h* w (;y))(z) = Jo exp(-zw) E*[exp(-zH^]R 2 , 2 (y,V2D^)^ b (dw) , 
in the sense of measurable functions, where the factor 

i? 2i2 (y,v / 2^)=V[l {T *< 5} ^( X |x 2 - y |)(2^)] 

= VDQ*(T b *<6) G m (2Dz) 
is well-defined and finite on {Re (z) > 0}. Here again recall the Section A. 3 function G a . 

6.4 The achievement time Laplace transform: This section establishes the Laplace 
transform of the normalized achievement time Hq by recalling ideas of [13, Appendix]. The 
principal finding is that the thus encoded excursion theoretic aspect of the problem leads 
to higher transcendental function denominators in the Laplace transform as follows: 

E[ex.p(-zH£)} = l/*(v / 2~Di) , Re (z) > , 

with \1/ the Section A. 2 function. The significance of this result at this stage is that it 
provides the missing factor in the Laplace transforms of Sections 6.2 and 6.3 and also 
shows all Laplace transforms to be finite on {Re (z) > 0}. 

The equivalent form of this result to be proved is the following key relation: 

E[exp(- \z 2 H*)] V(zVd) = tf(0) = 1, 

for any complex z with Re (z 2 ) > 0. Taking up the discussion of Section 6.1, the key 
relation is based on the Azema martingale [i. This process is defined in terms of the 
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Brownian meander of Section 6.1 by way of the equality: 

W* = m u (l) fj, u , ue[0,oo); 

it furnishes a martingale with respect to F + , the progressive enlargement of the Brownian 
filtration F by the sign of W* whose time-w step is given by F + (u) = J?~g u Vor(sgn W*), for 
any u > 0. Taking stochastic exponentials of this defining equality using the independence 
results of Fact 6.1 at time 1 there, 

£[exp {zW*-\z 2 t) | ^-+(g t )] = exp ( - \zH) *( Zf ji t ) , 

for any real z > 0. For fixed such z, the right-hand side of the equality at time t = Hq is 
a martingale by appealing to forms of the optional stopping theorem. The expectation of 
this martingale stopped at Hq is equal to its time-0 expectation which is equal to ^(0), 
and hence equals 1. Use the independence of H and W*(Hq) to obtain the key relation. 

7. Identification of the Cases I and II inversion problems 

We obtain the results of Section 3 in three steps from those of Sections 5 and 6 proceeding 
by reduction to the situation addressed by Theorem 4.1. Summarizing Sections 5 and 6, 
this section provides the first of these steps and identifies the functional relations to be 
considered. Since these transcribe the effects of two sources of stochasticity in two principal 
situations it is two representations which thus result as follows. 

Lemma 7.1: For any reals u > and y, we have in Case I: 

*<-'>-^(^)(s>)- 

Lemma 7.2: For any reals u > and y and with k in {1,2}, we have in Case II: 

Here Lemma 7.1 summarizes the findings of Section 6.2 and Lemma 7.2 those of Section 6.3, 
after a change of variables in both cases. 

8. Identification of the auxiliary random variables 

On comparison with Theorem 4.1 the results of Section 7 suggest immediate candidates 
for the functions R governing the auxiliary random variables An and Br to be intro- 
duced there. As a second step of the inversion procedure for the Section 3 results this 
section therefore identifies these functions and verifies their pertinent properties. These 
functions, it might be worth recalling from Section 6, transcribe the effects of two sources 
of stochasticity and their interrelations. In what follows let y denote an arbitrary real. 

Lemma 8.1: The function Ri(y, \fz) on the right-hand complex half-plane {Re (z) > 0} 
is obtained as a Laplace transform and Ri(y, y/z) = 0(z~ 3 ^ 2 ) as z tends to oo there. 

Lemma 8.2: We have S/ p ~ 1 (Ri(y, y/z)) = \/D Pb*,p( y ) ■ 

Proving these two results together we have to look at the functions R(z) = Ri(y, \fz) = 
v^Dexp (b* z)Gp( y )(z) where b* = b/\/D < 0, and we address their inversion right away by 
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replicating the functions p a , c of Section 3.1 with a = b* and c = (3(y). First let (3(y) < 0. 
Setting a(y) = — (b*+f3(y)), then appeal to Remark A. 3 o obtain: 

using Section A.l for inverting the first factor and Proposition A.l for inverting the second 
one. In the case /3(y) > 0, on the other hand, 

now using Proposition A. 2 for inverting the second factor. The proof of Lemma 8.2 is 
complete. Fom this discussion the asymptotic behaviour of # near oo required is immediate 
if b* or /3(y) are not 0. Otherwise, this behaviour is inherited from the one of Ni by way 
of Appendix B, and the proof of Lemma 8.1 is complete as well. 

Lemma 8.3: The function R2,i(y, \/z) on {Re (z) > 0} is obtained as a Laplace trans- 
form and #2,1(2/, y/z) = 0(z~ 3 ^ 2 ) as z tends to oo there. 

Lemma 8.4: We have i? 7-1 (#2,1(2/, \/z)) = \fT5 Po,p( y ) ■ 

These are special cases of Lemma 8.1 and 8.2 respectively; from the definitions in Sec- 
tions 3.3 and 6.3 recall how the functions #2,1 arise from functions Ri with b* = 0. 

Lemma 8.5: The function #2,2(2/, y/z) on {Re (2) > 0} is obtained as a Laplace trans- 
form and #2,2(2/, \/z) = 0(z~ 1 ) as z tends to 00 there. 

Lemma 8.6: We have (#2,2(2/, \/z))(u) = pt,(u,y) for any real u > 0. 

Setting R(z) = #2,2(2/, V^) reca11 R ( z ) = v^D J R S/ ? (x\ x _ y][ /^){z)ip b ^{x) dx from Sec- 
tion 6.3 in establishing these two results. 

Proof of Lemma 8.5 With the Section A.l Laplace transforms J5?(xa) defined on the right- 
hand half-plane, the function # is defined there as well. To determine its asymptotics there, 
express # by a direct computation in terms of the function Erfc as follows: 

sfzR{z) = V26D(V^/2)J2 exp(fc+£2/v / 2i)Erfc ((l/v / 25)(5v / 2i+£2/)) • 

Using the leading term of the asymptotic expansion of Erfc on the right-hand half-plane 
(see [9, Section 2.2]) the right-hand side of this expression is checked to behave like a scalar 
multiple of l/y/z as z tends to 00 there. Hence R(z) = 0(z~ 1 ) : as desired. 
To exhibit # as a Laplace transform, and to thus complete the proof, it is tempting to argue 
that its Laplace inverse is obtained by inversion under the sign of the defining integral. This 
turns out to be correct proceeding in two steps as follows. First delete an e-neighbourhood 
of y from the domain of integration of the defining integral, and effect Laplace inversion 
by inversion under the integral sign. Identify the integrand thus obtained as a continuous 
function bounded by an integrable function independent of y, and let e shrink to 0. 
Proof of Lemma 8.6 The last two step argument in proving Lemma 8.5 extends to the 
situation of Lemma 8.6. Reminding the Section A.l function cp a , we here obtain 

(^r) (r) = Erfc il dx = pb ^ y) ' 

where the last equality results noting Erfc (z/y/2Dr) = 2N 0i dt(z). The proof of Lem- 
ma 8.6 is complete as well. 
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9. Proof of the excursion density results 

Summarizing the development up to now, this section establishes the three convolution 
sum representation results, Theorems 3.1, 3.2 and 3.3, for the normalized excursion law 
and its two distinguished summands respectively. Referring to Section 7, each of these 
functions is the inverse of an explicitly given Laplace transform, and the principal idea 
is to effect these inversions analytically by applying Theorem 4.1. Apart from collecting 
terms the task at this point therefore reduces to two things. First, to establish for the 
numerators R of the Laplace transforms of Section 7 the asymptotic behaviour required in 
Theorem 4.1. Second, to identify the Laplace inverses of these numerators on division by 
the complex square root. This we address in turn in the three sections to follow. 

9.1 Proof of Theorem 3.1: To establish the Theorem 3.1 description of the nor- 
malized excursion density hi in Case I, we start from the Laplace inversion problem of 
Lemma 7.1. The idea is to apply Theorem 4.1 with the functions R(z) = Ri(y,^) = 
\/~D exp (b* z)Gp(y}{z) on {Re (z) > 0}. With the assumptions for this result satisfied by 
Lemma 8.1 indeed, 

K (u, y) = — ^ n=i 1 ( 3 >00) [—) 3T \- fT ) * !/„_! (— - - 

where the summation is over the finitely many positive integers n satisfying n < u/D. The 
Laplace inverses here are taken care of by Lemma 8.2 to be equal to Pb*,p( y ), and the proof 
of Theorem 3.1 is complete. 

9.2 Proof of Theorem 3.2: To establish the Theorem 3.2 description of the summand 
hi 2 l of the excursion law in Case II, we start from the inversion problem of Lemma 7.2 
with k = 1 there. Recall that this asks to integrate a Laplace inverse with respect to a 
density, and therefore we first look at the inversion problem alone. The idea is to apply for 
this inversion Theorem 4.1 with the function R(z) = i?2,i(y 5 V^)- With the assumptions 
of this result satisfied by Lemma 8.5 indeed, 

^ ,/ R(z) \, . v^oo (-1)™" 1 ifR(z)\ ( n\ 

\^(y/z)J y ' ^n=i (27r) n /2 l 2 >°°n > \ ) l \ 2/' 

for any r > 0, where only summands to indices n < 2r are not 0. The Laplace inverses 
here are then taken care of by Lemma 8.6 to be equal to Po,p( y ), as desired. 

A technical problem occurs on substitution of these results in Lemma 7.2. This is because 
one has to integrate there with respect to fit,(dw), and the integration variable enters into 
the respective number of the summands. Noting that these numbers are smaller than u/D, 
however, the representation of Theorem 3.2 is checked to follow, and the proof of this result 
is complete. 

9.3 Proof of Theorem 3.3: Establishing the Theorem 3.3 description of the summand 
^622 °f t ne excursion law in Case II essentially reduces to the problem considered in 
Section 9.1. In fact, writing out the second Laplace inversion problem of Lemma 7.2 with 
k = 2 there gives the representation: 
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where m u = min{5, u}. The Laplace inverses here are those of Sections 4.1 and 9.1 when 
b* = there, as formalized in Lemma 8.3 and 8.4. On inspection, substitution of the latter 
result's Laplace inverses yields the desired representation on harmonizing the number of 
summands as in Section 9.2. The proof of Theorem 3.3 is complete. 

10. Vista 

The paper has introduced a Parisian-style excursion law and determined its structure in 
the Brownian case, thus taking up and extending a development initiated in [13]. The 
characterization of our Brownian excursion law in terms of sums of independent random 
variables, however, does not lend itself readily to actual work with the law as it would 
be desirable for addressing the declared motivation from finance of [13], for example. 
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Figure 10.1. Comparison of Deltas of Down-and-In call [-] and Parisian Down-and-In call [o] 

Two methods for explicit handling of laws of sums of independent random variables are 
thus developed in the companion [12] to the present paper. Anticipating here some of the 
results there, for the Parisian barrier options proposed in [13] the methods are found to 
furnish effective and stable ways not just for valuation of the but also for hedging; this 
in particular so in situations where standard barrier options build up large Deltas 'near 
their barriers'; a typical example is featured in Figure 10.1 where Parisian barrier options 
permit a reduction of Deltas by some factor 5. We have thus come full circle, with the 
results of the present paper instrumental for this. 

Appendix A. Laplace transform pairs 

This appendix collects pertinent Laplace transform pairs. Here the Laplace transform is the 
linear operator 3 on the continuous functions of exponential type on (0, oo), the positive 
reals, defined as follows: it associates with any such function / the function 3(f) given by: 

&(f)( z ) = Jo°° exp (-zu)/(u) du , 
for any complex z in a half-plane contained sufficiently deep within the right-hand complex 
half-plane {z|Re(z) > 0}. The maps 3 ? (f) are analytic on such half-planes, and the 
operator 3 is an injection with inverse 3 , the inverse Laplace transform; see [3] or [4] 
for more detail. We moreover work with the principal branch of the complex logarithm on 
C \ (— oo, 0], the complex plane C cut along the non-positive reals (— oo, 0]. 

A.l The following three standard Laplace transforms on {z|Re (z) > 0} from [4, Beispiel 
8, p. 50f] originate with the heat equation: 



3(^a)(z) = exp(— otyfz) where "4> a (u) = 
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&(Xa)(z)= 6XP( ^ ) where Xa (u) = -}= exp ( - |-) , 

%)W = ^3 where y, a («) = Erfc (-^=) , 

for any real w > 0. Here Erfc(£) = {2/^/tt) /[^ )00 ) exp (— cc 2 ) c£r, for any complex £, is 
the complementary error function, and a is any complex with | arg(a)| < tv/4 such that 
Re (a) > for ip a . Hence the law ^ of the first passage time of Brownian motion to the 
level b is given on the real line R by: fib(dw) = ip by /2(w) dw. 

A. 2 As the first of two sets of functions to be considered define the functions N n on 
C \ (— oo, 0] for any integer n > by: 

N n = where N\(z) = V(-y/z)/Vz, 

for any z in C \ (— oo,0]. Following [16], the function \1> here is the generalization of the 
normal distribution given by the integral: 

ty(w) = J °° a; exp ( — ^x 2 +wx) dx , wEC. 

The second set is furnished by the functions v n = i>*^ on [0, oo) already used in Section 3. 
Let uq = be the Dirac delta function at and define v n for any integer n > 1 as an 
n-fold convolution on [0, oo) by: 

u n = v*( n) where v(u) = (2/Vtt) y/u/(2u+l) , u G [0, oo). 
Proposition A.l: We have Sfiy^) = N n on {Re (z) > 0}, for any integer n > 0. 

A. 3 This section concentrates on generalizations of the Section A. 2 function For any 
real a these are the functions F a and G a on C \ (— oo, 0] given by: 

F a (z) = Jq 00 xexp ( — \x 2 — \a — x\ y/z) dx , 

and 

G a {z) = F a {z)/yft. 

In terms of the Section A. 2 function h and the Section A.l complementary error function 
Erfc moreover define the function g a on [0, oo) by: 

9 Ju) = exp (- exp (- j^^) + ^Erfc 

for any u > 0. Then, in particular all G a are Laplace transforms as follows. 
Proposition A. 2: We have 2f{g a ) = G a on {Re (z) > 0}. 
Remark A. 3: If a < we have F a (z) = exp (ay/z)^(—y/z). 

Appendix B. Further properties of the function \1> 

This appendix develops pertinent properties of the function \1/ which from Section A. 2 for 
any complex number w is given by the integral: ^(w) = J °° xexp ( — |x 2 +w) dx. 

Developing the linear exponential factor of the integrand of \& in its series and integrating 
the resulting series term by term, yields the following series expansion: 

y( w ) = J2°°_ o anW n where a n = (2 n / 2 /n\) r(±(n + 2)). 
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This series is absolutely convergent for any complex number w, and its convergence is 
uniform on compact sets. As a first appplication it yields the key identity: 

\f/(w) = if>(-w) + y/2n w exp {\w 2 ) 

which connects the values of \1/ on the right-hand half-plane with those on the left-hand 
half-plane and has been noted in an equivalent form in establishing [1, Proposition 1 
point 2), p. 94]; the identity can also be obtained by partial integration of the defining 
integrals for ty(w) and ^>(—w) respectively by way of the partial integration identity: 

V2w) = 1 — y/wwexp (ty 2 )Erfc (w). 

This identity is the basis of the leading term expansion for \& on the left hand half plane: 

$ (-w) = 1/w 2 +R 2 (w) where \R 2 (w)\ < 6/\w\ 4 , 

for any complex w with Re (w) > 0, which is a special case of a general uniform asymptotic 
expansion of \1/ on the left-hand half-plane. 
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